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ABSTRACT 

Mechanism-based  models  for  the  evolution  of  defects  during  the 
thermomechanical  processing  of  aerospace  titanium-  and  nickel-base  alloys  are 
described.  These  defects  include  those  comprising  microstructural/metal-flow 
irregularities  and  those  that  are  damage-related  (i.e.,  cracks  and  cavities).  The 
development  of  undesirable/non-uniform  microstructures  and  cavities  during  the 
mill  processing  of  alpha/beta  titanium  alloys  is  addressed  first.  Relatively  simple, 
diffusion-based  models  of  spheroidization  and  coarsening  are  applied  to  quantify 
the  propensity  for  microstructure  non-uniformities.  Similarly,  first-order 
micromechanical  models  have  been  formulated  to  estimate  the  effect  of  local 
crystallographic  texture  on  non-uniform  flow,  the  generation  of  triaxial  stresses, 
and  cavity  growth/closure  in  alpha/beta  titanium  alloys  with  a  colony-alpha 
microstructure.  The  occurrence  of  non-uniform  grain  structures  (and  so-called 
ALA,  or  ‘as  large  as’,  grains)  in  cast,  wrought,  and  powder-metallurgy  superalloys 
are  also  discussed.  A  physics-based  model  to  treat  the  topology  of 
recrystallization  and  the  evolution  of  ALA  grains  in  such  materials  is  proposed. 
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1  Introduction 

The  thermomechanical  processing  of  aerospace  titanium-  and  nickel-base 
alloys  poses  a  number  of  challenges.  Often  intended  for  high-temperature 
applications,  these  materials  have  relatively  narrow  processing  windows  in  terms 
of  temperature  and  strain  rate  and  are  usually  produced  via  complex,  multi-step 
operations.  For  two-phase  alpha-beta  titanium  alloys  [1],  for  example,  ingots  with 
large  columnar  grains  are  subjected  to  hot  working  and  heat  treatment  in  the 
high-temperature  beta  (bcc)  phase  field  to  produce  a  finer  recrystallized  equiaxed 
grain  structure.  The  beta  phase  then  decomposes  during  cooling  below  the 
transus  temperature  (at  which  beta  -»  alpha  +  beta)  to  form  colonies  of  alpha 
platelets  (slow  cooling  rates)  or  a  basketweave  microstructure  (faster  rates) 
within  the  matrix  beta  grains.  Subsequent  hot  working  and  annealing  high  in  the 
two-phase  alpha/beta  field  is  used  to  breakup  and  spheroidize  the  transformed 
lamellar/basketweave  microstructure  and  thus  obtain  mill  product  with  a  structure 
comprising  fine  equiaxed-alpha  particles  in  a  matrix  of  beta. 

The  thermomechanical  processing  of  nickel-based  superalloys  also 
consists  of  multiple  steps  [2],  For  these  materials,  ingots  having  large  columnar 
grains  are  first  subjected  to  a  series  of  upsetting,  drawing,  and  reheating 
operations  above  the  solvus  temperature  (at  which  second  phases  such  as 
gamma  prime  are  dissolved)  in  order  to  obtain  a  recrystallized  equiaxed  gamma- 
grain  microstructure.  Secondary  hot  work  and  heat  treatments  below  the  solvus 
are  then  conducted  to  refine  the  grain  size  and  control  the  size  and  distribution  of 
second  phases. 
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Defects  that  are  developed  in  titanium-  and  nickel-base  alloys  are  most 
readily  grouped  according  to  whether  their  origin  is  related  to  melting,  primary 
processing,  secondary  (part)  processing,  or  final  heat  treatment.  Melt  processing 
of  titanium  ingots  can  lead  to  macrosegregation  of  oxygen,  beta  stabilizing 
elements,  etc.  that  produce  undesirable  variations  in  beta  transus  temperature 
and  hence  working  temperature  when  selected  relative  to  the  transus  [3],  In 
addition,  nitrogen-stabilized  hard  alpha  particles  and  high-density  inclusions  may 
be  retained  in  ingots  and  lead  to  substantial  property  debits  in  final  products  [4], 
Conversion  of  alpha/beta  titanium  ingots  to  billet  can  give  rise  to  microstructural, 
textural,  and  damage-type  defects.  These  include  alpha  lamellae  which  are  not 
fully  spheroidized  (i.e.,  retained  “spaghetti”  alpha),  regions  of  spheroidized  alpha 
particles  of  similar  texture  that  have  evolved  from  a  given  prior  colony  (leading  to 
“microtexture”),  macroscopic  defects  (e.g.,  surface  cracks  due  to  die  chill, 
centerbursts  due  to  the  development  of  secondary  tensile  stresses  and  poor 
workability),  and  microscopic  cavities  that  have  originated  at  the  prior-beta  grain 
boundaries  during  initial  alpha/beta  hot  working  [5,  6], 

The  secondary  processing  of  titanium  parts  can  also  result  in  a  variety  of 
defects,  many  of  which  are  metal  flow  in  nature  and  driven  by  improper  die 
design  and/or  material  properties  [5],  These  include  laps  (metal  folding  over 
itself),  flow-through  defects  (e.g.,  metal  flow  past  a  rib  in  a  closed  die  without 
filling  it),  pipe  (e.g.,  suck-in  defects  at  the  end  of  extrudates),  and  shear  bands. 
Shear  bands  can  be  especially  troublesome  because  they  may  originate  due  to 
numerous  factors  such  as  poor  die/preform  design,  die  chilling,  friction,  improper 
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process  parameters  (e.g.,  high  strain  rates  that  give  rise  to  substantial 
deformation  heating),  and  material  characteristics  (e.g.,  large  degrees  of  flow 
softening).  Part  processing  may  also  give  rise  to  defect  generation  during  final 
heat  treatment.  For  instance,  the  development  of  non-uniform  texture  throughout 
a  titanium  part  during  deformation  processing  can  produce  non-uniform  beta 
grain  size  during  final  beta  annealing  due  to  the  interaction  of  grain  growth  and 
texture  [7], 

A  number  of  the  defects  found  in  wrought  titanium  alloys  also  can  occur  in 
wrought  nickel-base  superalloys.  However,  there  are  also  some  which  are 
specific  to  superalloys  such  as  microsegregation,  ALA  (‘as  large  as’)  grains, 
abnormal  grains,  and  thermal  cracking  [8],  Microsegration  results  from  alloying 
element  partitioning  during  solidification  and  the  sluggish  kinetics  of  the 
homogenization  of  a  number  of  alloying  elements  utilized  in  superalloys  to 
optimize  service  properties.  Grain-size  non-uniformities  may  result  from  partial 
recrystallization  (ALA  grains)  or  the  rapid,  preferential  growth  of  a  small  subset  of 
so-called  ‘abnormal’  grains  during  supersolvus  heat  treatment.  Last,  gross 
cracking  may  occur  during  heating  which  is  too  rapid  (e.g.,  ingot/billet  reheating 
prior  to  hot  working)  or  quenching  of  thick-section  parts  following  final  solution 
treatment.  In  both  cases,  the  development  of  temperatures  gradients/thermal 
stresses  act  in  concert  with  limited  ductility  to  cause  fracture  [9], 

The  great  breadth  of  defects  that  can  be  developed  in  aerospace  alloys 
precludes  an  extensive  treatment  of  each  in  the  present  paper.  Thus,  the 
objective  of  the  present  work  is  to  provide  an  overview  of  several  relatively 
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straightforward,  physics-based  approaches  for  analyzing  defect  generation  and 
illustrate  the  application  of  each  with  an  actual  production-type  problem.  The 
methods  to  be  discussed  comprise  the  bulk-diffusion  analysis  of  the  evolution  of 
microstructural  defects  in  titanium  alloys,  plasticity  analyses  for  cavitation  during 
hot  working  of  textured  titanium  alloys,  and  a  mesoscale  model  for  describing  the 
non-uniformity  of  recrystallization  during  the  hot  working  and  heat  treatment  of 
superalloys. 

2  Microstructural  Defects  in  Alpha/Beta  Titanium  Alloys 

The  breakdown  of  colony  (or  basketweave)  alpha  to  obtain  a  fine, 

equiaxed-alpha  microstructure  (Figure  1)f  during  primary  processing  of 
alpha/beta  titanium  alloys  in  the  two-phase  phase  field  is  controlled  by  a  number 
metallurgical  processes.  Each  of  these  processes  has  an  effect  on  the  final  alpha 
particle  size  and  morphology  and  thus  on  final  strength  and  fatigue  resistance, 
among  other  properties.  First,  preheating  prior  to  hot  working  leads  to  coarsening 
of  the  alpha  platelets;  i.e.,  their  thickness  increases  and  length  decreases  (Figure 
2a).  Spheroidization  during  subsequent  deformation  occurs  by  shearing-type  or 
boundary  splitting  mechanisms  [10].  For  this  reason,  prior  static  coarsening  of 
the  alpha  platelets  results  in  substantial  and  undesirable  increases  in  the  strains 
required  for  dynamic  spheroidization  [1 1],  These  strains  are  typically  in  excess  of 
those  that  can  be  imposed  during  industrial  processes  used  to  fabricate  mill 
products  in  sizes  required  for  secondary  processes  such  as  part  forging.  Hence, 
the  completion  of  spheroidization  must  be  accomplished  during  static,  post- 

t  In  backscattered-electron  (BSE)  images  taken  in  a  scanning-electron  microscope  (SEM),  the 
alpha  phase  is  dark,  and  the  beta  (or  martensitic  alpha)  phase  is  white  (or  gray). 
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deformation  heat  treatment  (Figure  2b).  Static  spheroidization  occurs  by 
boundary  splitting  and,  to  a  larger  degree,  on  termination  migration  [12,  13]. 
Finally,  if  the  time  for  completion  of  the  majority  of  the  static  spheroidization  is 
exceeded,  measurable  static  coarsening  of  the  equiaxed  alpha  particles  can 
occur  (Figure  2c).  As  with  other  metallic  materials,  static  coarsening  leads  to  the 
growth  of  large  alpha  particles  and  the  dissolution  of  small  particles. 

2.1  Diffusion  Analyses  of  Microstructure  Evolution.  Each  of  the 
various  static  heat-treatment  processes  that  control  microstructure  evolution 
during  the  breakdown  of  the  colony-alpha  morphology  in  alpha/beta  titanium 
alloys  (i.e. ,  coarsening  of  alpha  platelets,  spheroidization  via  termination 
migration,  and  coarsening  of  a  spheroidized  microstructure)  is  driven  by  a 
reduction  in  overall  alpha/beta  interface  energy.  Furthermore,  it  has  been  found 
that  the  controlling  mechanism  by  which  surface  area/energy  is  reduced  is  the 
bulk  diffusion  of  solutes  through  the  beta  matrix  [14,  15].  In  particular,  solutes 
diffuse  from  regions  of  high  chemical  potential/concentration  to  regions  of  lower 
chemical  potential/concentration.  For  the  simplest  case  in  which  the  two  phases 
are  ideal,  terminal  solid  solutions,  the  Gibbs-Thompson  equation  gives  the  solute 
concentration  of  the  matrix  (Cp(R))  next  to  a  precipitate  (referred  to  as  a)  with 
principal  radii  of  curvature  equal  to  Ri  and  R2,  viz., 

Cp(R)/Cp  (R=oo)  =  1  +  (YapQ/kBT)  ( J_  +  J_ ) ,  (1 ) 

Ri  R2 
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in  which,  Cp(R  =  qo)  denotes  the  composition  next  to  a  flat  precipitate  (zero 
curvature),  yaP  is  the  alpha/beta  interface  energy,  Q  is  the  atomic  volume,  kB  is 
Boltzmann’s  constant,  and  T  is  the  absolute  temperature. 

The  application  of  Equation  (1)  is  simplest  for  the  spheroidization  of  a 
lamellar  platelet  (assumed  to  have  a  pancake  shape)  via  termination  migration 
[14].  In  such  an  instance,  solute  is  transferred  from  the  periphery  to  the  flat  (plan) 
surfaces  of  each  (assumed  non-interacting)  precipitate.  As  suggested  by 
Courtney  and  Malzahn  Kampe  for  a  similar  2D  problem  [16],  a  relatively  simple 
application  of  Fick’s  Law  to  solute  diffusion  at  the  beginning  and  intermediate 
points  of  the  process  can  provide  reasonable  estimates  of  the  spheroidization 
time.  For  example,  in  the  beginning  of  the  diffusional  process  during  which  the 
plan  surfaces  are  still  relatively  flat,  the  concentration  difference  between  the 
edges  and  plan  surfaces  is  as  follows: 

[Cp  (R=oo)](yapQ/kBT)  ( 2  +  — 1 — )  .  (2) 

t  w  t  w 

2~+4 

Here,  t,  w  are  the  platelet  thickness  and  diameter,  respectively.  Geometrical 
considerations  provide  estimates  of  the  diffusion  distance  (and  hence 
concentration  gradient)  and  the  area  through  which  the  solute  flux  passes. 

Semiatin,  et  at.  [14}  extended  the  Courtney  and  Malzahn  Kampe 
spheroidization  analysis  to  a  pancake  geometry  and  modified  it  to  account  for  the 
fact  that  the  alpha  and  beta  phases  in  titanium  alloys  are  not  terminal  solid 
solutions.  The  final  result  for  the  spheroidization  time,  TVd,  was  the  following: 
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^vd 

X' 


(3a) 


n3-[o.328n7/3(1+A/1_o.763n-4/3)2] 

2(1 + n)  |  o.5n1/3+o.665n2/3  ^ 
3(0.5-0.572r|-1/3)  3(0.143 +  0.934n“1/3) 


in  which 


r\  =  (w/t)  +  0.5 

t  =  t3RT/DpCFyapVM 


CF  =  Composition  Factor  = 


Cp(l-Cp) 

(Ca-Cp)2[1  +  dlnr/dlnCp] 


(3b) 

(3c) 

(3d) 


In  Equations  (3),  R  is  the  gas  constant,  Dp  is  the  diffusivity  of  the  rate  limiting 
solute  in  beta  titanium,  Vm  is  the  molar  volume,  and  Cp  (=Cp(R  =  oo))  and  Ca  are 
the  equilibrium  solute  concentrations  (expressed  as  atomic  fractions)  of  the  rate- 
limiting  solute  in  the  beta  and  alpha  phases,  respectively. 

The  other  two  diffusion  problems  mentioned  above  are  more  complex.  For 
platelet  coarsening,  solute  is  transferred  from  the  edges  of  a  given  platelet  to  its 
own  plan  surfaces  or  to  those  of  adjacent  precipitates.  The  added  geometric 
complexity  relative  to  that  for  the  analysis  of  spheroidization  of  an  individual 
platelet  is  presently  under  investigation.  For  static  coarsening  of  a  distribution  of 
spherical  precipitates,  a  modified  Lifshitz-Slyosov-Wagner  (MLSW)  analysis  [17, 
18]  has  been  found  to  describe  kinetics  for  alpha/beta  titanium  alloys  such  as  Ti- 
6AI-4V  [15], 


2.2  Validation  of  Static  Spheroidization  Model.  The  validity  of  Equation 
(3)  to  quantify  the  spheroidization  behavior  and  persistence  of  remnant  lamellae 
was  established  using  observations  for  Ti-6AI-4V  heat  treated  at  955°C  (Figure 
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2b).  At  this  temperature,  the  shapes  of  the  lamellae  that  existed  after  the 
completion  of  boundary  splitting  (time  «  1h)  approximated  the  idealized  pancake 
geometry  assumed  in  the  diffusion  model.  For  Ti-6AI-4V,  the  diffusion  of 
vanadium  through  the  beta  matrix  is  rate-limiting  (in  comparison  to  the  diffusion 
of  aluminum). 

Model  predictions  of  tVC|/t’  and  hence  xvd  for  spheroidization  of  the 

remnant  lamellae  are  given  in  Table  I.  The  predicted  values  for  the 
spheroidization  time  for  the  pancake-shape  geometry  ranged  from  approximately 
5  to  15  hours  for  5  of  the  7  platelets;  the  other  two  lamellae  were  predicted  to 
require  times  of  the  order  of  30  hours.  These  values  of  spheroidization  time 
compare  well  with  the  measured  time  of  12-14  hours  to  achieve  approximately  90 
volume  pet.  of  spheroidized  microstructure,  as  defined  by  alpha  particles  with  an 
aspect  ratio  of  less  than  2:1  [12,  14].  Thus,  it  is  not  surprising  that  a  small  volume 
fraction  of  modest  aspect-ratio,  partially-spheroidized  lamellae  would  remain  after 
14  hours  (Figure  2b). 

3  Plasticity  Analyses  for  Cavitation  in  Alpha/Beta  Titanium  Alloys 

Internal  cavitation  represents  a  second  major  class  of  defects  in 

alpha/beta  titanium  alloys.  The  control  of  cavitation  is  very  important  because 
such  damage  may  lead  to  poor  service  properties  as  well  as  premature  failure. 

Cavitation  during  hot  working  and  superplastic  forming  usually  occurs  via 
nucleation,  growth,  and  coalescence  processes.  Nucleation  represents  perhaps 
the  most  difficult  problem  to  model,  and,  depending  on  the  mechanism  by  which 
it  occurs,  different  approaches  have  been  proposed.  For  particle-containing 
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materials,  such  as  aluminum  and  nickel-base  alloys,  several  techniques  have 
been  formulated  to  predict  the  critical  particle  size  above  which  cavity  nucleation 
is  likely  to  occur  at  a  given  strain  rate  due  to  the  inability  to  relax  local  stress 
concentrations  via  diffusion  [19-21],  In  titanium  alloys,  which  typically  do  not 
contain  second-phase  particles,  cavitation  is  most  common  during  the 
breakdown  hot  working  of  the  colony-alpha  microstructure  (Figure  3).  In  such 
cases,  an  alternate  analysis,  assuming  nucleation  due  to  slip  intersections  with 
and/or  vacancy  coalescence  at  the  prior-beta  grain  boundaries  followed  by  initial 
cavity  growth  under  constrained-plasticity  conditions,  has  been  proposed  [22], 

Following  nucleation,  cavity  growth  usually  occurs  via  a  general-plasticity- 
controlled  mechanism.  For  uniaxial-tension  conditions,  growth  of  an  isolated 
cavity  is  then  usually  described  by  the  following  relation  [23]: 

dr  it  ,  .x 

dWr  ■  (4) 

in  which  r  denotes  the  particle  radius  as  a  function  of  imposed  strain  s,  and  r\  is 
the  cavity-growth  parameter. 

In  the  subsections  that  follow,  several  analyses  for  the  growth  ( and 
shrinkage)  of  cavities  in  a  typical  alpha/beta  titanium  alloy  (Ti-6AI-4V)  with  a 
colony-alpha  microstructure  are  summarized. 

3.1  Analysis  for  the  Average  Cavity  Size.  The  sizes  of  the  cavities  that 
develop  during  hot  working  are  frequently  affected  by  coalescence  with  adjacent 
cavities  as  well  as  continuous-nucleation  phenomena  that  occur  as  deformation 
proceeds.  To  quantify  these  phenomena,  Nicolaou  and  Semiatin  [24,  25] 
conducted  a  number  of  mesoscale  simulations  of  cavitation  during  uniaxial 
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tension  and  found  that  the  dependence  of  the  average  cavity  radius  r  on  true 
axial  strain  s  could  be  fit  by  an  expression  of  the  following  form: 

£  =  „r(H  +  Cv)  (5) 

in  which  Cv  denotes  the  cavity  volume  fraction,  and  H  is  a  constant  that  depends 
on  the  cavity-nucleation  conditions.  For  the  case  of  an  initial  cavity  array  for 
which  growth  and  coalescence  occur  with  no  further  nucleation,  it  was  found  that 
H  =  1/3.  On  the  other  hand,  if  cavity  nucleation  occurs  continuously  in  addition  to 
growth  and  coalescence,  the  value  of  H  was  shown  to  decrease  to  0.2. 

Two  modifications  to  Equation  (5)  are  required  to  gage  the  evolution  of  the 
average  cavity  size  under  complex  stress  conditions  [26],  First,  the  axial  strain  s 
is  replaced  by  the  effective  strain  s .  Second,  the  cavity-growth  parameter  for 
uniaxial  tension  (rj)  is  replaced  by  a  cavity-growth  parameter  r|ts  which  is  a 
function  of  stress  triaxiality,  or  the  ratio  of  the  mean  stress  ctm  to  the  effective 
stress  a.  Using  hot  torsion  and  notched  tension  tests,  Nicolaou,  et  al.  [27-29] 
determined  the  ratio  of  r\ls/r\.  These  results  were  bounded  by  those  from  the 
theoretical  analysis  of  Rice  and  Tracey  [30]  for  cavity  growth  under  cold-working 
conditions  and  the  semi-empirical  relation  derived  by  Pilling  and  Ridley  [31]  for 
superplastic  forming  (Figure  4). 

The  approach  developed  by  Nicolaou  and  Semiatin  [26]  to  predict  the 
average  cavity  size  developed  under  complex  stress  conditions  was  validated 
using  conventional  hot  pancake  forging  of  Ti-6AI-4V.  During  such  a  forging 
operation,  barreling  gives  rise  to  secondary  tensile  stresses  (and  thus  cavitation) 
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whose  magnitude  decreases  with  distance  from  the  free  surface  (Figure  5a). 
Thus,  continuum  finite-element-method  (FEM)  predictions  of  the  local  stress 
state  (Figure  5a),  a  critical  damage  factor  for  cavity  nucleation,  and  values  of  r|ts 
were  combined  with  Equation  (5)  to  predict  the  variation  of  the  average  cavity 
size  with  distance  from  the  free  surface.  The  predictions  showed  good 
agreement  with  measurements  (Figure  5b). 

3.2  Micromechanical  Analysis  for  the  Size  of  the  Largest  Cavities. 

The  classical  plasticity-controlled  cavity-growth  analysis  is  incapable  of  predicting 
the  size  of  the  largest  cavities  that  are  developed  during  hot  working  of  highly- 
anisotropic  materials  such  as  titanium  alloys.  In  such  cases,  the  dependence  of 
flow  properties  on  crystallographic  orientation  must  be  taken  into  account.  To  this 
end,  Bieler,  et  al.  [32]  and  Nicolaou,  et  al.  [26]  developed  an  approximate 
micromechanical  analysis  for  the  stress  state  and  strain  partitioning  between 
adjacent  “hard”  and  “soft”  colonies  in  alpha/beta  titanium  alloys,  i.e.,  those 
colonies  with  either  large  or  small  Taylor  factors,  Mh  or  Ms.  One  of  the  basic 
assumptions  of  the  model  was  that  the  cavity  formed  at  the  boundary  between 
the  hard  and  soft  colonies  was  small  relative  to  the  grain  size  and  hence  had  little 
influence  on  the  stresses  and  strains  that  would  be  developed  in  the  absence  of 
the  cavity. 

For  the  case  of  a  complex  macroscopic  stress  state,  such  as  would  be 
imposed  during  metalworking,  the  micromechanical  analysis  comprised  yield 
criteria  for  the  hard  and  soft  colonies,  equilibrium  conditions,  and  a  self- 
consistent  estimate  of  the  strain  partitioning  between  the  hard  and  soft  colonies. 


12 


The  yield  criteria  in  terms  of  the  axial  (az),  radial  (ar),  and  hoop  (a0)  stress 
components  were  the  following: 


ah  =_Tf{(ae  _GZ  )2+K  _ar  )2+(ar  _a0  )2^/2  (6a) 

V2  h  h  h  h  h  h 

a  =— jL{(a  -a  )2+(a  -a  )2+(a  -a  )2}^2  (6b) 

sA/olv0  z  7  v  z  r  7  v  r  0  7  J 
^  s  s  s  s  s  s 

in  which  the  subscripts  h  and  s  refer  to  the  hard  and  soft  colonies,  and  asand  ah 

denote  the  flow  stress  of  the  soft  and  hard  colonies,  respectively. 

Cavities  were  assumed  to  open  on  colony/grain  boundaries  perpendicular 
to  the  hoop  direction.  In  such  cases,  the  macroscopic  axial  (a^ac)  and  radial 

(o™ac )  stresses  determined  from  FEM  analysis  were  each  assumed  to  be  a  rule- 

of-mixtures  average  of  the  corresponding  stress  components  in  the  hard  and  soft 
colonies,  thus  satisfying  the  equilibrium  considerations,  i.e., 

azac  =fhazh  +(1_fh)az  (7a) 

^aC=fh\+(1-fhHs  •  (7b) 

in  which  fh  denotes  the  volume  fraction  of  hard  colonies.  Furthermore,  it  was 
assumed  that  the  hoop  stress  was  the  same  in  the  soft  and  hard  colonies  and 
was  equal  to  the  macroscopic  (FEM)  stress  a™0 ,  i.e., 


0.mac  =  (j  =  (j 
ue  °eh 


(7c) 


If  the  flow  stresses  of  the  hard  (a.)  and  the  soft  (a  )  colonies  are  known, 


the  axial  and  radial  stresses  (and  hence  stress  triaxiality)  in  each  colony  can  be 
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determined  by  solving  Equations  (6)  and  (7).  Because  flow  stress  is  rate 
sensitive  at  hot  working  temperatures,  as,  ah,  and  the  corresponding  strain  rates 

were  determined  by  applying  a  self-consistent  model  [33],  In  this  approach,  the 
flow  stress  of  the  aggregate  amac  is  given  by  the  following  equation: 

amac=fhah+(1-fh)as  (8) 

Assuming  that  the  flow  stress  is  strain  rate  and  strain  sensitive,  Equation  (8)  can 
then  be  rewritten  as 

amac  =fhkh(th)m(8h)n  +  (1-fh)ks(ts)m(8s)n  (9) 

in  which  kh,  ks  are  the  strength  coefficients  of  the  hard/soft  colonies,  which 
depend  only  on  the  Taylor  factors  Mh/Ms;  s,  s  denote  the  effective  strain  and 
strain  rate,  respectively;  and  m  and  n  are  the  strain  rate  sensitivity  index  and  the 
strain  exponent.  Using  the  self-consistent  model,  the  values  of  strain  rates  (and 
therefore  strains)  in  the  hard  and  soft  colonies  can  then  be  determined  for 
specific  values  of  kh  and  ks. 

For  the  same  hot  pancake  forging  example  discussed  in  Section  3.1,  the 
micromechanical  model  was  used  to  predict  the  stress  triaxiality  and  the  size  of 
the  cavities  between  adjacent  hard  and  soft  colonies,  assuming  fh  =  0.6  (Figure 
6).  Model  predictions  of  the  size  of  the  largest  cavities  showed  good  agreement 
with  observations  for  Mh/Ms  ~  2,  or  values  typical  of  the  actual  texture  [32], 

Several  of  the  important  aspects  of  the  micromechanical  model  for  the 
size  of  the  largest  cavities  developed  during  hot  working  of  titanium  have  been 
validated  by  detailed  crystal-plasticity  FEM  (CPFEM)  calculations.  For  example, 
CPFEM  [34]  has  shown  that  the  strain  partitioning  between  hard  and  soft 
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colonies  is  well  replicated  well  by  the  self-consistent  method  of  Semiatin,  et  al. 
[33],  Similarly,  the  stress  triaxiality  in  the  soft  colony  estimated  by  the  simple 
equilibrium  approach  shows  reasonable  agreement  with  CPFEM  results.  For 
example,  CPFEM  has  indicated  that  the  triaxiality  in  the  softer  and  harder 
colonies  (near  the  hard-soft  colony  interface)  is  approximately  1.0  or  0.05, 
respectively,  for  Mh/Ms  =  2.5,  in  reasonable  agreement  with  the  equilibrium 
calculations  (Figure  6a). 

3.3  Cavity-Shrinkage  Analysis.  Cavities  developed  during  the  hot 
working  of  anisotropic  materials  such  as  alpha/beta  titanium  alloys  may  be 
decreased  in  size  or  closed  by  changes  in  strain  path  which  give  rise  to  stress 
triaxiality  that  is  compressive  in  nature.  For  instance,  cavity  shrinkage  for  Ti-6AI- 
4V  subjected  to  forward  torsion  followed  by  reversed  torsion  or  uniaxial 
compression  is  illustrated  in  Figure  7. 

An  analysis  to  describe  cavity-shrinkage  kinetics  has  been  recently 
developed  and  validated  by  Nicolaou  and  Semiatin  [35,  36],  In  their  approach, 
the  stresses  developed  within  the  soft  colony/grain  (during  the  strain-path  change 
which  gives  rise  to  compressive  triaxiality)  are  estimated  in  the  same  manner  as 
in  the  cavity-growth  analysis  described  in  Section  3.2.  To  estimate  the  strain-rate 
components  and  thus  the  densification  rate  (sn  +  s22  +  s33 ),  the  stresses  are  then 
inserted  into  a  hybrid  model  for  the  consolidation  of  porous  media  [37],  i.e., 

KW  £  V(]+v)a:.  +(i  -2v)am8ij]  (10) 
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in  which  sy  denotes  the  strain  rate  tensor,  <jy  is  the  deviatoric  stress  tensor,  ctm 

the  mean  stress,  a  is  the  effective  stress,  v  is  the  Poisson's  ratio  of  the  porous 
body;  and  (j>  is  the  stress  intensification  factor.  The  function  K(D)  is  associated 
with  the  relative  density  D;  it  is  commonly  assumed  to  be  equal  to  D  for  relative 
densities  greater  than  90  pet.  of  full  density. 

Cavity-shrinkage  predictions  from  the  analysis  Nicolaou  and  Semiatin  [35, 
36]  showed  good  agreement  with  observations  (e.g.,  the  results  in  Figure  8  for 
compression  following  torsion).  In  particular,  the  analysis  correctly  quantified  the 
more  rapid  closure  kinetics  during  compression  compared  to  reversed  torsion 
(each  following  forward  torsion).  Such  a  result  was  explained  on  the  basis  of  the 
higher  levels  of  compressive  stress  triaxiality  in  compression  and  the  orientation 
of  cavities  in  torsion  relative  to  the  applied  stresses,  among  other  factors. 

4  Mesoscale  Model  for  Recrystallization  of  Superalloys 

Grain  size/shape  irregularities  comprise  another  major  class  of  defects 

found  in  aerospace  alloys  and  are  of  special  concern  for  iron-  and  nickel-base 
superalloys.  The  prediction  of  such  defects  requires  models  that  treat 
deformation,  recrystallization,  and  grain-growth  during  thermomechanical 
processing.  A  number  of  techniques  have  been  developed  to  describe  these 
phenomena.  For  example,  the  cellular-automata  (CA)  approach  can  be  used  to 
simulate  microstructure  evolution  on  a  relatively  small  scale,  but  cannot  be 
applied  at  every  node  of  an  FEM  mesh  for  a  finite-size  billet/workpiece  due  to 
computation  limitations.  By  contrast,  simple,  phenomenological  techniques  (such 
as  that  based  on  the  Avrami  formulation)  do  not  have  this  drawback,  but  typically 
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lack  the  sophistication  required  to  treat  complex  problems  involving  strain  rate 
and  temperature  transients.  Hence,  an  intermediate  (meso-)  scale  model  for 
microstructure  evolution  that  is  mechanism  based  (as  in  the  CA  technique)  and 
computationally  efficient  (and  thus  applicable  to  every  node  of  an  FEM  mesh) 
has  been  developed  recently  [38-40]  and  is  currently  being  applied  to  analyze  the 
occurrence  of  grain-structure  defects  in  superalloys. 

4.1  Formulation  of  Mesoscale  Model.  The  mesoscale  modeling 
technique  comprises  two  main  parts,  a  geometrical  framework  and  one 
incorporating  driving-force  relations  that  describe  the  actual  evolution  of 
microstructure.  To  minimize  computation  requirements,  the  geometric  framework 
does  not  represent  every  grain  in  an  aggregate.  Rather,  the  microstructure 
description  relies  on  averaging  over  grain  populations  that  are  characterized  by 
similar  properties  and  which  are  expected  to  exhibit  similar  evolution.  This 
requirement  leads  to  the  definition  of  grain  families  referred  to  as  meso-structure 
units,  or  MSUs.  To  deal  with  recrystallization,  at  least  two  MSUs  are  necessary: 
one  for  the  initial  grains  and  another  for  the  recrystallized  ones.  In  the  case  of 
multi-hit  deformation,  such  as  during  cogging  of  ingots,  at  least  two  generations 
of  recrystallized  grains  are  required:  one  for  the  grains  that  appeared  during 
previous  hits,  which  may  have  grown  and  then  stored  some  energy  again,  and 
those  that  nucleated  during  the  current  or  latest  deformation  step. 

The  geometric  description  of  microstructure  also  relies  on  several 
assumptions  in  terms  of  grain  shape  (e.g.,  spherical  vs  ellipsoidal)  and  the 
location  of  nucleation  sites  (near  existing  grain  boundaries  or  inside  the  initial 
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grains,  as  in  particle-stimulated  nucleation  (PSN)  of  recrystallization).  Geometric 
relations  deduced  from  these  assumptions  (primarily  statistical  expectancies  of 
the  surface  of  contact  between  grains)  constrain  the  geometric  variables  of 
MSUs  (e.g.,  grain  size,  grain  density)  into  a  framework  that  evolves  in  response 
to  two  kinds  of  inputs:  nucleation  rates  and  grain-boundary  velocities.  From  these 
inputs,  the  geometric  framework  allows  the  autonomous  and  internally-consistent 
calculation  of  the  geometric  evolution  of  the  various  grain  populations  during  their 
growth -related  interactions.  A  few  variables  are  added  to  each  MSU  to  quantify 
dislocation  density,  sub-boundary  density/misorientation,  or  intrinsic  properties 
such  as  its  Taylor  factor. 

Driving-force  equations  are  used  to  evaluate  the  evolution  of  the  MSU 
parameters  and  to  deduce  from  them  the  required  inputs  of  the  geometric 
framework.  The  most  important  of  these  equations  describes  the  dislocation 
mean  free  path  (strain  hardening),  dynamic  and  metadynamic  recovery,  sub¬ 
boundary  generation  rate,  the  kinetics  of  sub-boundary  disorientation,  nuclei  size, 
grain  boundary  mobility  and  velocity,  etc.  Details  on  the  mesoscale  model  can  be 
found  in  References  39  and  40. 

4.2  Mesoscale  Model  Predictions  for  ALA  Grains.  The  mesoscale 
model  has  been  applied  to  treat  the  evolution  of  ALA  grains  during  the 
breakdown  of  superalloy  ingots  with  coarse,  columnar  starting  grains  [40],  Such 
materials  usually  have  a  grain  size  of  ~5  mm  diameter  x  ~50  mm  length  and 
contain  -1000  PSN  sites  per  mm3  which  are  clustered  in  groups  of  -10  particles 
[41]. 
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Two  modeling  approaches  have  been  used  for  the  ALA  problem,  one  in 
which  the  discrete  zones  of  grain-boundary  (“necklace”)  and  intragranular  (PSN- 
recrystallized)  grains  persist  throughout  the  simulation  and  the  other  in  which  the 
discrete  necklace  and  PSN  zones  are  converted  to  a  solely  wrought-necklace 
microstructure  when  the  volume  fraction  of  intragranular  PSN  zones  reaches  2/3 
(Figure  9a).  It  has  been  found  that  the  evolution  of  the  dynamically-recrystallized 
volume  fraction  as  a  function  of  strain  is  essentially  identical  for  the  two  cases. 
However,  the  benefit  of  the  second  approach  is  evident  when  examining  the 
predicted  size  of  the  remnant  unrecrystallized  (ALA)  regions  (Figure  9b).  For  the 
first  case,  the  initial  ingot  grains  seem  to  disappear  suddenly.  Hence,  this 
approach  in  general  fails  to  provide  physical  insight  into  the  evolution  of  ALA 
grains.  On  the  other  hand,  useful  information  about  the  size  of  the  remnants  of 
the  initial  grains  is  obtained  from  the  second  modeling  approach  involving  the 
topological  conversion.  At  the  moment  when  the  initial  ingot  grains  are  converted 
into  wrought-like  grains,  the  curve  describing  the  size  of  the  remnants  exhibits  a 
change  because  the  information  it  provides  is  no  longer  based  on  the  initial  ingot 
grains  but  rather  the  volume  of  the  new  wrought-like  remnants.  The  boundaries 
of  these  new  grains  are  inherited  primarily  from  the  interfaces  of  the  percolating 
intragranular  zones  in  which  most  of  the  recrystallization  has  occurred.  Thus, 
they  are  the  true  remnants  of  the  initial  ingot  grains.  Following  the  conversion, 
the  curve  describing  the  size  of  the  remnant  unrecrystallized  areas  exhibits  a 
smooth  decrease  until  recrystallization  is  complete  (Figure  9b). 
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5  Future  Challenges 

Physics-based  models  are  invaluable  for  the  design  and  control  of 
thermomechanical  processes  for  aerospace  alloys.  However,  a  number  of  major 
challenges  remain  for  the  full  application  of  such  models  to  minimize  or  eliminate 
defects  from  finished  components.  These  include  the  following: 

•  Micromechanical  understanding  of  phenomena  such  as  cavity  nucleation 
during  hot  working  and  abnormal  grain  growth  during  final  beta  heat  treatment 
of  titanium  alloys. 

•  Refined  CPFEM  models  which  treat  the  generation,  multiplication,  and 
storage  of  dislocations.  Such  models  are  needed  for  the  construction  of 
related  physics-based  models  of  dynamic,  metadynamic,  and  static 
recrystallization  of  metallic  alloys  and  the  defects  associated  with  such 
metallurgical  phenomena. 

•  Models  for  residual  stress  development  and  distortion,  including  internal- 
state-variable  constitutive  models  for  complex  thermal  histories. 

•  Input  data  for  models,  including  thermodynamic/kinetic  data  in  both  the  solid 
and  liquid  states  (e.g.,  phase  equilibria,  mobility/diffusivity  data),  grain¬ 
boundary  properties  (energy,  mobility),  and  large-strain,  high-strain-rate 
constitutive  data  and  models  at  cold,  warm,  and  hot-working  temperatures. 
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Table  I.  Geometry  of  Remnant  Alpha  Lamellae  and  Model  Predictions  of 
Spheroidization  Time  at  for  Ti-6AI-4V  at  955°C 


Lamella 

ID 

Lamella 
Diameter, 
w  (pm) 

Lamella 
Thickness,  t 
(pm) 

h 

Tvd/'t 

Tvd  (h) 

A 

16.4 

2.7 

6.5 

1.39 

6.1 

B 

17.7 

3.2 

6.1 

1.13 

7.8 

C 

18.2 

3.6 

5.5 

0.84 

8.6 

D 

21.8 

3.6 

6.5 

1.39 

14.4 

E 

24.1 

2.3 

11.1 

6.16 

15.6 

F 

27.3 

3.6 

8.0 

2.52 

26.1 

G 

33.2 

2.3 

15.1 

13.6 

34.4 
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Figure  Captions 
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Figure  5. 


Microstructures  developed  in  Ti-6AI-4V  during  primary  processing: 
(a)  optical  and  (inset)  SEM  BSE  micrographs  of  the  colony-alpha 
microstructure  and  (b)  SEM  BSE  micrograph  of  the  fine,  equiaxed- 
alpha  microstructure. 

Diffusion-controlled  processes  during  the  breakdown  of  the  colony- 
alpha  microstructure:  (a)  static  coarsening  of  alpha  plates,  (b)  static 
spheroidization  of  remnant  alpha  plates  following  hot  working,  and  (c) 
classical  static  coarsening  of  equiaxed-alpha  particles. 

Cavitation  developed  during  hot  tension  testing  of  Ti-6AI-4V  with  a 
colony-alpha  microstructure:  (a)  optical  micrograph,  (b)  high- 
magnification  SEM  BSE  image,  and  (c)  electron-backscatter- 
diffraction  (EBSD)  inverse-pole  figure  map  indicating  presence  of 
hard  (basal-oriented)  and  soft  colonies  around  cavity.  The  tension 
axis  is  vertical  in  all  micrographs. 

Ratio  of  the  cavity-growth  parameter  under  non-uniaxial  stress  states 
to  that  for  uniaxial  tension  (r|ts/r|)  as  a  function  of  the  stress  triaxiality 
(ratio  of  the  mean  stress  cjm  to  the  effective  stress  a ).  The  data 
points  (determined  using  hot  torsion  and  notched  tension  tests  on  Ti- 
6AI-4V)  [27-29]  are  compared  to  predictions  from  models  of  Rice  and 
Tracey  [30]  and  Pilling  and  Ridley  [31]. 

Model  results  for  cavitation  in  a  Ti-6AI-4V  pancake  forging:  (a) 
continuum  FEM  model  predictions  of  the  local  stress  state  at  the  mid- 
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height  as  a  function  of  effective  strain  and  distance  from  the  free 
surface  and  (b)  comparison  of  measurements  and  predictions  of  the 


average  cavity  size  after  a  height  reduction  of  50  pet.  at  815°C.  In  (b), 
the  sensitivity  of  the  model  predictions  to  input  data  is  illustrated  by 
varying  r|ts  by  ±  25  pet  relative  to  its  nominal  value  (i.e. ,  Q  =  0.75,  1, 
or  1.25). 

Figure  6.  Micromechanical-model  predictions  for  cavitation  in  Ti-6AI-4V  with  a 
colony-alpha  microstructure  which  was  hot  pancake  forged  at  815°C: 
(a)  stress  triaxiality  in  adjacent  hard  and  soft  colonies  at  the 
equatorial  free  surface  as  a  function  of  effective  strain  and  the  Taylor- 
factor  ratio  (Mh/Ms)  and  (b)  cavity  size  after  a  50-pct.  reduction  as  a 
function  of  the  Taylor-factor  ratio  and  distance  from  the  equatorial 
free  surface.  The  cavity-size  predictions  in  (b)  are  compared  to 
measurements  (data  points). 

Figure  7.  Micrographs  illustrating  the  effect  of  strain  path  on  cavitation  in 
samples  of  Ti-6AI-4V  with  a  colony-alpha  microstructure  subjected  to 
hot  deformation  at  815°C  and  an  effective  strain  rate  of  0.04  s'1:  (a) 
forward  torsion  to  a  surface  effective  strain  of  0.99,  or  forward  torsion 
to  a  surface  effective  strain  of  0.99  followed  by  (b)  reversed  torsion  to 
a  surface  effective  strain  of  0.25,  or  (c)  uniaxial  compression  along 
the  prior  torsion  axis  to  a  height  strain  of  -0.29.  The 
torsion/compression  axis  is  vertical  and  the  radial  direction  is 
horizontal  in  the  micrographs. 
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Figure  8.  Model  predictions  of  the  reduction  in  cavity  volume  fraction  as  a 
function  of  macroscopic  effective  strain  during  compression  of  Ti-6AI- 
4 V  with  a  colony-alpha  microstructure  following  torsion  to  the 
indicated  surface  effective  strains.  Torsional  prestraining  and 
subsequent  compression  were  both  conducted  at  815°C  and  an 
effective  strain  rate  of  0.04  s'1.  The  predictions  are  compared  to 
experimental  results  (data  points). 

Figure  9.  Mesoscale  model  for  simulating  the  evolution  of  recrystallization 
nucleated  at  grain  boundaries  and  intragranular  particles  in 
superalloy  ingots  with  coarse,  columnar  grains:  (a)  schematic 
illustration  of  the  topology  of  recrystallization  and  (b)  predicted  size  of 
the  remnant,  un recrystallized  (ALA)  grains  for  the  cases  in  which  (i) 
the  original  topology  was  retained  throughout  the  simulation  (solid 
line)  and  (ii)  the  original  topology  was  converted  to  a  necklace-only 
one  once  the  intragranular  recrystallized  regions  had  percolated 
throughout  the  initial  ingot  grains  (broken  line).  The  nucleation  rate 
was  1  nucleus  per  100  pm2  of  boundary  and  per  unit  strain,  the  grain¬ 
boundary  velocity  was  10  pm  per  unit  strain,  and  the  volume  of  each 
nucleus  was  1000  pm3. 
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Figure  1.  Microstructures  developed  in  Ti-6AI-4V  during  primary  processing: 

(a)  optical  and  (inset)  SEM  BSE  micrographs  of  the  colony-alpha 
microstructure  and  (b)  SEM  BSE  micrograph  of  the  fine,  equiaxed- 
alpha  microstructure. 
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(b) 


Figure  2.  Diffusion-controlled  processes  during  the  breakdown  of  the  colony- 
alpha  microstructure:  (a)  static  coarsening  of  alpha  plates,  (b)  static 
spheroidization  of  remnant  alpha  plates  following  hot  working,  and  (c) 
classical  static  coarsening  of  equiaxed  alpha  particles. 
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Figure  3.  Cavitation  developed  during  hot  tension  testing  of  Ti-6AI-4V  with  a 
colony-alpha  microstructure:  (a)  optical  micrograph,  (b)  high- 
magnification  SEM  BSE  image,  and  (c)  electron-backscatter- 
diffraction  (EBSD)  inverse-pole  figure  map  indicating  presence  of 
hard  (basal-oriented)  and  soft  colonies  around  cavity.  The  tension 
axis  is  vertical  in  all  micrographs. 
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Figure  4.  Ratio  of  the  cavity  growth  parameter  under  non-uniaxial  stress  states 
to  that  for  uniaxial  tension  (r|t7r|)  as  a  function  of  the  stress  triaxiality 
(ratio  of  the  mean  stress  cjm  to  the  effective  stress  a ).  The  data 
points  (determined  using  hot  torsion  and  notched  tension  tests  on  Ti- 
6AI-4V)  [27-29]  are  compared  to  predictions  from  models  of  Rice  and 
Tracey  [30]  and  Pilling  and  Ridley  [31]. 
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Figure  5.  Model  results  for  cavitation  in  a  Ti-6AI-4V  pancake  forging:  (a) 
continuum  FEM  model  predictions  of  the  local  stress  state  at  the  mid¬ 
height  as  a  function  of  effective  strain  and  distance  from  the  free 
surface  and  (b)  comparison  of  measurements  and  predictions  of  the 
average  cavity  size  after  a  height  reduction  of  50  pet.  at  815°C.  In 
(b),  the  sensitivity  of  the  model  predictions  to  input  data  is  illustrated 
by  varying  r|ts  by  ±  25  pet  relative  to  its  nominal  value  (i.e.,  Q  =  0.75, 
1,  or  1.25). 
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Figure  6.  Micromechanical-model  predictions  for  cavitation  in  Ti-6AI-4V  with  a 
colony-alpha  microstructure  which  was  hot  pancake  forged  at  815°C: 
(a)  stress  triaxiality  in  adjacent  hard  and  soft  colonies  at  the 
equatorial  free  surface  as  a  function  of  effective  strain  and  the 
Taylor-factor  ratio  (Mh/Ms)  and  (b)  cavity  size  after  a  50-pct.  reduction 
as  a  function  of  the  Taylor-factor  ratio  and  distance  from  the 
equatorial  free  surface.  The  cavity-size  predictions  in  (b)  are 
compared  to  measurements  (data  points). 
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Figure  7.  Micrographs  illustrating  the  effect  of  strain  path  on  cavitation  in 
samples  of  Ti-6AI-4V  with  a  colony-alpha  microstructure  subjected  to 
hot  deformation  at  815°C  and  an  effective  strain  rate  of  0.04  s'1:  (a) 
forward  torsion  to  a  surface  effective  strain  of  0.99,  or  forward  torsion 
to  a  surface  effective  strain  of  0.99  followed  by  (b)  reversed  torsion  to 
a  surface  effective  strain  of  0.25,  or  (c)  uniaxial  compression  along 
the  prior  torsion  axis  to  a  height  strain  of  -0.29.  The 

torsion/compression  axis  is  vertical  and  the  radial  direction  is 
horizontal  in  the  micrographs. 


Cavity  Fraction  (Pet) 


Figure  8.  Model  predictions  of  the  reduction  in  cavity  volume  fraction  as  a 
function  of  macroscopic  effective  strain  during  compression  of  Ti-6AI- 
4V  with  a  colony-alpha  microstructure  following  torsion  to  the 
indicated  surface  effective  strains.  Torsional  prestraining  and 
subsequent  compression  were  both  conducted  at  815°C  and  an 
effective  strain  rate  of  0.04  s'1.  The  predictions  are  compared  to 
experimental  results  (data  points). 
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Figure  9.  Mesoscale  model  for  simulating  the  evolution  of  recrystallization 
nucleated  at  grain  boundaries  and  intragranular  particles  in 
superalloy  ingots  with  coarse,  columnar  grains:  (a)  schematic 
illustration  of  the  topology  of  recrystallization  and  (b)  predicted  size  of 
the  remnant,  un recrystallized  (ALA)  grains  for  the  cases  in  which  (i) 
the  original  topology  was  retained  throughout  the  simulation  (solid 
line)  and  (ii)  the  original  topology  was  converted  to  a  necklace-only 
one  once  the  intragranular  recrystallized  regions  had  percolated 
throughout  the  initial  ingot  grains  (broken  line).  The  nucleation  rate 
was  1  nucleus  per  100  pm2  of  boundary  and  per  unit  strain,  the  grain¬ 
boundary  velocity  was  10  pm  per  unit  strain,  and  the  volume  of  each 
nucleus  was  1000  pm3 
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